Cosmological Constraints on an Invisibly Decaying Higgs 
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Working in the context of a proposal for collisional dark matter, we derive bounds on the Higgs 
boson coupling g' to a stable light scalar particle, which we refer to as phion {(p), required to solve 
problems with small scale structure formation which arise in coUisionless dark matter models. We 
discuss the behaviour of the phion in the early universe for different ranges of its mass. We find 
that a phion in the mass range of 100 MeV is excluded and that a phion in the mass range of 1 GeV 
requires a large coupling constant, g' > 2, and < 130 GeV in order to avoid overabundance, in 
which case the invisible decay mode of the Higgs boson would be dominant. 



PACS numbers: 95.35.+d, 98.62., 14.80.B 



I. INTRODUCTION 

We have recently discussed the role of a Ught, sta- 
ble, strongly self-coupled scalar field in solving the prob- 
lems of the Cold Dark Matter (CDM) model for struc- 
ture formation in the Universe, concerning galactic scales 
[01 . Our proposal involves a particle physics- motivated 
model, where the DM particles are allowed to self-interact 
so as to have a large scattering cross section and negli- 
gible annihilation or dissipation. The self-interaction re- 
sults in a characteristic length scale given by the mean 
free path of the particle in the halo. This idea was origi- 
nally proposed to suppress small scale power in the stan- 
dard CDM model and has been recently revived, 
in a general context, in order to address the issues dis- 
cussed above 0|. Our model [|] is a concrete realization 
of this idea, which involves an extra gauge singlet as the 
self-interacting, non-dissipative cold dark matter parti- 
cle. Following Ref. ||], we call this scalar particle phion, 
(j), and assume that it couples to the Standard Model 
(SM) Higgs boson, h, with a Lagrangian density given 
by: 

C = lid,^r~lmlcj^^-^^^^+g'vcl^^h , (1) 

where g is the phion self-coupling constant, its mass, 
V = 246 GeV is the Higgs vacuum expectation value and 
g' is the coupling between (p and h. A model along these 
lines have been previously discussed Clearly the in- 
teraction term between the phion and the Higgs boson 
arises from a quartic interaction ^(jP'H'^ , where H is the 
electroweak Higgs doublet. As shown in 0, the (j) mass 
does not arise from spontaneous symmetry breaking since 



this would yield a tiny scalar self-coupling constant. The 
phion mass in should be regarded phenomcno- 
logical parameter arising from a more encompassing the- 
ory. 

As is well known, scalar particles have been repeatedly 
invoked as DM candidates [0-0; however, our proposal 
has the salient feature that it brings about a connection 
with the SM Higgs boson which could arise in extensions 
of the SM. For instance, the hidden sector of heterotic 
string theories does give rise to astrophysically interest- 
ing self-interacting scalars For reasonable values of 
g' , the new scalar would introduce a novel invisible decay 
mode for the Higgs boson. This could, in principle, pro- 
vide an explanation for the failure in finding the Higgs 
boson at accelerators sofar which will be tested at 
future colliders. 

On the astrophysical front, recent observational data 
on large scale structure, cosmic microwave background 
anisotropics and type la supernovae suggest that ^tot ~ 
1, of which ^baryons ~ 0.05 and r^A « 0.65 jl^; the 
remaining contribution, ^dm ~ 0.3 (apart from neu- 
trinos that may contribute a small fraction) , comes from 
dark matter (DM) , which determines the hierarchy of the 
structure formation in the Universe. The most promi- 
nent theories of structure formation are now ACDM and 
QCDM, which consist, respectively, of the standard Cold 
Dark Matter (CDM) model supplemented by a cosmo- 
logical constant or a dark energy, i.e. a negative pressure 
component. 

In the CDM model, initial Gaussian density fluctu- 
ations, mostly in non-relativistic coUisionless particles, 
the so-called cold dark matter, grow during the inflation- 
ary period of the Universe and evolve, via gravitational 
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instability, into the structures one observes at present. 
However, it has been found that the CDM model cannot 
sucessfully accomodate the data observed on all scales. 
For instance, N-body simulations predict a number of 
halos which is a factor ^ 10 larger than the observed 
number at the level of Local Group |l6yi7| ]. Further- 
more, astrophysical systems that are DM dominated, e.g. 
dwarf galaxies |lq-|2Ct| and low surface brightness galax- 
ies show shallow matter-density profiles with finite 
central densities. This is in contradiction with galactic 
and galaxy cluster halos in high resolution N-body simu- 
lations ]2^-|25t, which have singular cores, with p ~ r~'^ 
and 7 in the range between 1 and 2. This can be under- 
stood as cold collisionless DM particles do not have any 
characteristic length scale leading, due to hierarchical 
gravitational collapse, to dense dark matter halos with 
negligible core radius. 

It is relevant to stress that recent numerical simula- 
tions [p6|-p9| indicate that self-interaction of DM parti- 
cles does bring noticeable improvements on properties of 
the CDM model on small scales. We point out, however, 
that a numerical simulation that takes into account the 
salient features of our proposal is still missing. 

At present, (p particles are non-relativistic, with typical 
velocities v ~ 200 km s^^, and, therefore, it is impossible 
to dissipate energy creating more particles in reactions 
like (fxf) —^ (j)(f)(j)(j3. Thus, as only the elastic channel is 
kinematically accessible, near threshold, the cross section 
is given by: 
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16tts 647rm3, 



(2) 



A limit on and g can be obtained by demanding 
that the mean free path of the particle 0, A^, is in the 
interval 1 kpc < < 1 Mpc B. Hence, we have: 
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(3) 



where and are, respectively, the number and mass 
density of </> particles in the halo. Eqs. and (||) imply: 

= 2.1 X 10^ f ^) 

VGeV/ \MpcJ 



0.4 GeVcm" 



GeV~ 



(4) 



which, in turn, leads to: 
= 13 g2/3 ( 
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VMpcy \0.4GeVcm' 



MeV 



(5) 



In what follows, we shall see how the requirement that 
il^/i^ ~ 0.3, i.e. that the phion is a suitable DM can- 
didate, and that the phion is able to explain small scale 
structure, leads to bounds on the couphngs g and g'. 



II. PHION DENSITY ESTIMATE 

In the ensuing discussion, we shall bear in mind the 
most recent lower bound on the mass of the SM Higgs 
boson, as it emerges from the combined LEP data gath- 
ered at energies up to 205.9 GeV, mh > 112.3 GeV at 
95% confidence level On the other hand, SM preci- 
sion data indicates that m/j < 188 GeV at 95% confidence 
level [|o) and nih < 306 GeV at 99% confidence level 

If g' is sufficiently small, phions decouple early in the 
thermal history of the Universe and are diluted by sub- 
sequent entropy production. In Ref. j^], we have con- 
sidered out-of-equilibrium phion production via inflaton 
decay in the context of A'^ = 1 Supergravity inflationary 
models (see. e.g. |3^ and references therein). 

On the other hand, for certain values of the coupling g', 
it is possible that </> particles are in thermal equilibrium 
with ordinary matter. In order to determine whether this 
is the case, we will make the usual comparison between 
the thermalization rate Tth and the expansion rate of the 
Universe H. 

The thermalization rate is given by 



Tth = n < (JannVrel > 



(6) 



where n = 1.2 x T^/tt'^ is the density of rclativistic phions 
and < (TannVrei > IS the annihilation cross section aver- 
aged over relative velocities. On the other hand, the 
expansion rate is given by: 
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47r3g, 
45 
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1.66 X gl^^ 



rp2 



(7) 



At temperatures above the electroweak phase transi- 
tion. Tew — 300 GeV, a typical value in many extensions 
of the SM where one hopes to find the required features to 
achieve successful baryogenesis [^,0 , the order param- 
eter (the vacuum expectation value of the Higgs field) 
vanishes, v(T) = 0, and therefore the (j)(j)h coupling is 
non-operative. However, as mentioned above, this inter- 
action term has its origin in the 4-point coupling, cfxjjhh, 
which can bring, at high temperatures, the phion-Higgs 
system into thermal equilibrium. Using the temperature 
as the center-of-mass energy, the cross section is given 
by: 
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which implies that phions are in thermal equilibrium for 
temperatures smaller than 
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Therefore, phions will never be in thermal equilibrium 
before the electroweak phase transition if g' < 10~^. 

Thermal equilibrium can be achieved just below Tewi 
when the trilinear coupling is operative, if g' is such that 
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the thermalization rate, Tth, is larger than the Hubble 
expansion rate. We now compute the required bounds 
on g'. 

The phion annihilation cross section at high energies 
(T > rrih) is given by a relativistic Breit-Wigner reso- 
nance formula: 



l2^2 , 



(10) 



where Th is the total Higgs decay rate. At the resonance 
peak (s = TO^) it becomes 



47r 

iVrel = -^BRQl 



Using 



(11) 



(12) 



we obtain the decoupling temperature in the limit mh ^ 



150 



(13) 



This implies that, in order to have a decoupling tempera- 
ture of the order of the Higgs mass, the coupling constant 
should be fairly small: 



g' ~ 10 



-10 



(14) 



where we have introduced the SM value oiVh = 3.2 MeV, 
obtained from the code HDECAY |^^, used to compute 
the width of a 115 GeV Higgs boson. 

Therefore, if g' > 10"^", the phions will be brought 
into thermal equilibrium right after the electroweak 
phase transition. If this is the case, there are two pos- 
sible scenarios depending whether they decouple while 
relativistic or otherwise. 

In order to study these scenarios, we need to determine 
the decoupling temperature at T ~ <C Tew, in which 
case the phion annihilation cross section involves virtual 
Higgs exchange (h*), as in Figure |l|, and is given by [ p6[ : 
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ann ^re. 



Vrel 



{4ml " "^ft)^ + "^l^l 



F 



X 



where 



Fx = lim 



(15) 
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and Th*x denotes the width for the decay h* ^ X (X 

since we are dealing with inelastic scattering only), 
for ruh' = 2m^. For the mass range of interest to us, 
TO^ ~ 10 - 100 MeV, one has Fx ~ 10"" 

In this case, the relationship between the coupling con- 
stant g' and the decoupling temperature Tjy is given by 




g,b,W-,Z,... 



FIG. 1. Feynman diagramm for phion annihilation via 
Higgs exchange. 



,2 ^ 5 J. (m JlOOGeV)^ 



(Td/McV) 



(17) 



If ff' ^ 0.1, the phions decouple while relativistic and 
are as abundant as photons. Since we are interested in 
stable light phions, it is a major concern avoiding phion 
overproduction if it decouples while relativistic. In fact, 
in this case, there is an analogue of Lee- Weinberg limit 
for neutrinos (see e.g. pH]): 



0.08 



(18) 



which results in a very stringent bound, < 4 eV, 
and implies that the phion self-coupling constant should 
satisfy g < 2.5 x 10~^°, in order to solve the small scale 
structure problem that exists with conventional collision- 
less cold dark matter. 

In order that the phions decouple non-relativistically 
and their abundance reduces to acceptable levels with- 
out fine-tuning the self-coupling constant, g' > 0.1 is re- 
quired. In this case, using standard methods to compute 
the phion relic abundance ||3q,p8|, one obtains: 



1.07 X W^XF 
gj Mp{ 



(19) 



where g^ denotes the number of degrees of freedom in 
equilibrium at annihilation and xp = m^/Tp is the in- 
verse of the freeze-out temperature in units of the phion 
mass. The relevant crosss section is the phion annihila- 
tion cross section invol ving virtual Higgs exchange, Eq. 
(p!5|), with Fx ~ 10"^'^ | |37| |. The freeze-out temperature 
is determined by the solution of the equation 



Xp 



\n[Q.{m{g,xp)-^/^Mpm^{aannVrei)] ■ (20) 



In order to get xp > 1 and apply Eq. (|l|), one must 
have g' > 1.2, for = 50 MeV and rrih = 115 GeV. 
However, due to the smallness of the cross section, the 
relic abundance of the phion is several orders of mag- 
nitude larger than the observed value. Hence, one con- 
cludes that a phion in this mass range is excluded. 
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FIG. 2. Contour of n^h'^ = 0.3 as a function of m?i (in 
GeV) and g' , for = 0.5 GeV (top), 1.0, 1.5 and 2 GeV 
(bottom). 

The smallness of the phion annihilation cross section 
for ~ 50 MeV is due to the smah factor Fx — 10^^^. 
However, this factor increases significantly with larger 
phion masses. Considering ~ 1 GeV, then it follows 
that Fx — 10"''. We find that the requirement U^h^ ~ 
0.3 implies that > 500 MeV and 5' > 2, a solution 
which holds only for m^j < 130 GeV. Heavier phion 
and Higgs particles tend to make fi^/i^ > 0.3. We have 
verified that our solutions respect the condition xp > 3, 
meaning that the phion behaves as a cold dark matter 
candidate. Of course, our solutions are at the edge of 
validity ofperturbation theory. Our results are depicted 
in Figure g. 

For these large values of the coupling constant, the 
decay width of the Higgs into phions is given by: 

Tih -> cbcj^) = 5.23 ( YY^^) g'' GeV. (21) 

Hence, in this case the Higgs width is totally domi- 
nated by the invisible decay mode and this model can be 
easily tested at future colliders (see e.g. ]39|). 



as discussed in [y. In this case, the phion does not 
contribute to the invisible Higgs boson decay channel. 
For g' > 10^1", we have found that, if g' < 0.1, the 
phion decouples while still relativistic and a limit for its 
mass, 7710 < 4 eV can be derived, which, in turn, implies 
in a strong bound on the phion self-coupling constant, 
g < 10~^, if the phion is required to solve the CDM 
model problems on small scales. On the other hand, 
if 5' > 1, the phion decouples while non-relativistic; 
however, its abundance is not cosmologically acceptable 
for phion masses in the range of 50 — 100 MeV due to 
the small annihilation cross section. For masses in the 
range of 0.5 — 2 GeV, we have found that abundances of 
ri^/i^ ~ 0.3 require large values of the coupling g' ~ 2.5 
and nih < 130 GeV. In this scenario, the Higgs width 
is dominated by the invisible h — > (j)(f) mode and can be 
tested at future colliders. 

Regarding the origin of the phion mass, it is possible 
that it arises from the cancellation between a tachyonic 
mass, due to the dynamics of fields in a more encompass- 
ing theory together with the contribution from the quar- 
tic interaction, g''^v'^. In fact, in case the only relevant 
fields are the phion and the Higgs boson, the phion mass 
is dominated by a term which arises from electroweak 
symmetry breaking, namely, ~ g'v, and if the phion 
is to solve the small structure problem, this would imply 
g' ~ 5 X 10~^ which, in turn, brings us to the overpro- 
duction problem. 

Finally, we comment on the recent observation ^] 
that CDM problems at small scales can be fixed 
only through s-channel annihilation of DM particles, 

that is (JannVrel OC COUSt SO that (JannVrel = 2.5 X 

10~^(TO/GeV) GeV~^. We have verified that these con- 
ditions cannot be met in our model. On more general 
grounds, it is difficult to see how substantial amounts of 
annihilating DM could be obtained from any processes 
in the early Universe, given that the required cross sec- 
tion is so large, although it is clearly a quite interesting 
challenge to build concrete particle physics models that 
exhibit such features. 
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duced by out-of-equilibrium decay of the inflaton field. 
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